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This  research  program  utilizes  satellite  and  surfaced-derived  cloud  observations  together 
with  standard  meteorological  measurements  to  evaluate  and  improve  our  ability  to 
accurately  diagnose  cloud  coverage.  Results  of  this  research  will  be  used  to  compliment 
existing  or  future  parameterizations  of  cloud  effects  in  global  and  regional-scale 
meteorology  forecast  models,  since  nearly  all  cloud  parameterizations  must  specify  a 
fractional  area  of  cloud  coverage  when  calculating  radiative  or  dynamic  cloud  effects, 
and  current  parameterizations  rely  on  rather  crude  cloud  cover  estimates. 
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1.  Project  Statement  of  Work 

This  research  program  utilizes  satellite  and  surface-derived  cloud  observations  together 
with  standard  meteorological  measurements  to  evaluate  and  improve  our  ability  to 
accurately  diagnose  cloud  coverage.  Results  of  this  research  will  be  used  to  compliment 
existing  or  future  parameterizations  of  cloud  effects  in  global  and  regional-scale 
meteorology  forecast  models,  since  nearly  ail  cloud  parameterizations  must  specify  a 
fractional  area  of  cloud  coverage  when  calculating  radiative  or  dynamic  cloud  effects, 
and  current  parameterizations  rely  on  rather  crude  cloud  cover  estimates. 

During  the  first  phase  of  this  research  program,  our  goal  is  to  evaluate  and  improve  the 
methods  for  calculating  cloud  cover  within  a  mesoscale  meteorology  model.  To 
accomplish  this,  a  mesoscale  meteorology  model  will  be  quantitatively  evaluated  using 
the  U.  S.  Air  Force  3DNEPH  and  RTNEPH  satellite-derived  cloud  fields.  Hourly- 
averaged  distributions  of  the  model-derived  cloud  fields  will  be  compared  with  observed 
clouds  at  the  finest  spatial  and  temporal  resolution  of  the  corresponding  datasets. 
Algorithms  currently  used  in  global  or  mesoscale  meteorological  models  to  assess  cloud 
coverage  will  be  objectively  evaluated.  We  will  then  carry  out  an  innovative  search  for 
relationships  between  cloud  coverage  and  numerous  meteorological  factors  such  as 
relative  humidity,  stability,  wind  shear,  moisture  convergence,  precipitation  rate  and  other 
parameters. 


During  the  second  phase  of  this  research,  the  cloud  cover  data  and  improved 

parameterizations  of  cloud  coverage  developed  during  the  first  phase  will  be  incorporated  _ 
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into  a  mesoscale  meteorology  model.  Model  forecasts  which  utilize  the  observed  cloud 
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coverage  and  depth  should  be  improved  relative  to  forecasts  which  crudely  specify  cloud  •^mn-i^ir  t-< 
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2.  Project  Progress: 


1  Oct.  1991  -  30  Sept.  1992 


During  the  first  year  of  this  research  program,  we  have  compiled  and  reviewed  a  list  of 
formulations  used  by  various  research  groups  to  specify  cloud  cover.  We  find 
considerable  variability  between  formulations  used  by  various  climate  and  meteorology 
models,  and  under  some  conditions,  one  formulation  will  produce  a  zero  cloud  amount, 
while  an  alternate  formulation  calculates  95%  cloud  cover  under  the  same  environmental 
conditions.  All  formulations  hypothesize  that  cloud  cover  is  predominantly  determined 
by  the  average  relative  humidity,  although  some  formulations  allow  local  temperature 
lapse  rates  and  vertical  velocities  to  influence  cloud  amount.  All  formulations  specify 
that  cloud  fraction  is  zero  below  a  “threshold”  relative  humidity,  and  there  are  several 
estimates  of  this  vertically-varying  “threshold”  based  on  physical  arguments  or  “tuning” 
of  free  parameters  to  match  global  albedo  estimates.  We  find  that  a  majority  of  the 
formulations  are  based  on  limited  observations,  and  evaluations  of  algorithm  accuracy  are 
based  on  longer-term  compilations  of  cloud  cover  observations  that  cannot  resolve  the 
vertical  distribution  of  cloud  amount. 

In  order  to  explicitly  evaluate  the  accuracy  of  these  and  future  cloud  cover  algorithms,  we 
have  prepared  an  evaluation  database  of  cloud  cover  and  related  meteorological 
observations  during  a  springtime  midlatitude  cyclone.  This  growing  dataset  contains 
simultaneous  observations  of  cloud  cover  fraction,  moisture,  wind  speeds,  temperatures  at 
over  5000  (320  km)2  areas  over  the  northeast  U.  S.  at  15  tropospheric  levels.  Cloud 
cover  is  derived  from  the  U.  S.  Air  Force  3DNEPH  satellite  archive,  and  related 
meteorology  is  extracted  from  the  NMC  analysis.  We  have  developed  initial  versions  of 
a  set  of  analysis  programs  designed  to  search  for  and  quantify  correlations  between  cloud 
cover  and  other  meteorological  factors.  This  cloud  observational  database  is  being  used 
to  assess  the  accuracy  of  existing  cloud  cover  algorithms. 


Progress  Report:.  Walcsh 


0306  4 


We  find  that  most  cloud  cover  algorithms  tend  to  underpredict  cloud  coverage  at  any 
particular  level  in  the  troposphere,  especially  under  dry  conditions,  when  we  find  10  -  20 
%  cloud  cover  at  relative  humidities  as  low  as  20  -  30%  within  certain  tropospheric  layers 
averaged  over  (320  km)^  areas.  We  also  find  that  cloud  amount  appears  to  increase  with 
increasing  humidity  at  ail  relative  humidities,  in  stark  contrast  to  existing  formulations, 
which  contain  hq  cloud  cover  sensitivity  to  relative  humidity  at  humidities  below  the  70  - 
90  %  “critical”  humidities  used  by  these  models.  We  see  no  clear  evidence  that  cloud 
cover  vanishes  at  a  “critical  humidity”.  Based  on  these  initial  comparisons,  we  feel  that 
current  short-term  meteorology  models  cannot  adequately  assess  the  changes  in  cloud 
cover  that  may  results  from  small  changes  in  relative  humidity  induced  by  larger-scale 
forcings.  Additionally,  a  potentially  important  feedback  between  climate  change  and 
changes  in  cloud  cover  are  probably  not  adequately  simulated  by  current  models  of  global 
warming. 

We  have  discovered  an  apparent  discrepancy  between  observed  relative  humidity  in  the 
upper  troposphere  under  dry  conditions,  and  calculations  of  relative  humidity  using  the 
NCAR  Mesoscale  meteorology  model,  with  the  model  calculating  excessive  relative 
humidity  in  the  upper  troposphere  under  unstable  conditions.  This  is  most  likely 
attributable  to  an  incorrect  parameterization  of  moisture  redistribution  by  the  convective 
parameterization  used  by  the  model.  Such  an  error  would  lead  to  an  inaccurate  cloud 
coverage  calculation  if  it  were  present  in  a  forecast  model. 

We  have  begun  development  of  an  improved  convective-scale  vertical  mixing 
parameterization  for  use  in  our  mesoscale  meteorological  model  in  order  to  alleviate  the 
deficiency  in  the  NCAR  mesoscale  meteorological  model.  At  this  point,  the  model  is 
being  tested  as  a  stand-alone,  single  column  model  that  is  being  tested  with  observations 
of  tropical  convection.  This  convective  mixing  algorithm  will  be  based  on  a  detraining 
cloud-plume  description  of  buoyantly-driven  cloud  motions,  and  include  a  reasonably 
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sophisticated  description  of  multi-stream  cloud  up-  and  down-drafts  with  somewhat  more 
sophisticated  microphysics  than  are  currendy  incorporated  into  most  cumulus 
parameterizations.  This  new  cumulus  parameterization  will  be  used  to  assess  the  heating 
and  moistening  tendencies  due  to  buoyantly-induced  convection  within  hydrostatic 
models  of  the  atmosphere  that  cannot  explicitly  compute  these  effects. 

In  addition  to  the  time  devoted  by  Dr.  Walcek  to  this  research  effort,  two  post-doctoral 
research  associates  (Dr.  Kesu  Zhang  and  Dr.  Hu  Qi)  and  a  graduate  student  are 
contributing  to  various  technical  aspects  of  this  project.  We  have  also  purchased 
microcomputer  and  workstation-level  computer  software  and  hardware  as  needed  by  this 
project. 

Dr.  Walcek  presented  results  from  this  research  effort  at  two  scientific  meetings:  the  11th 
International  Conference  on  Clouds  and  Precipitation  held  in  Montreal,  Canada,  17-21 
August  1992;  and  the  WMO  Cloud  Microphysics  and  Applications  to  Global  Climate 
Change  Workshop  held  in  Toronto,  Canada  from  10-14  August  1992.  We  have  also 
prepared  a  manuscript  for  journal  publication  summarizing  the  research  results  described 
above.  Meeting  abstracts  citations  are  listed  at  the  end  of  this  progress  report,  and  a  copy 
of  the  reviewed  manuscript  (undergoing  minor  revisions)  is  provided  with  this  report. 

4.  Budget  Status: 

As  of  30  September  1992, 98%  of  the  moneys  budgeted  for  the  first  year  have  been 
expended  approximately  as  initially  budgeted.  Moneys  budgeted  for  computer  expenses 
were  used  to  purchase  microcomputer  and  workstation  hardware  and  software.  Results  of 
the  initial  analysis  of  cloud  cover  were  presented  at  international  meetings  in  Canada,  and 
travel  expenses  associated  with  these  meetings  were  included  with  the  original  budget. 
Costs  associated  with  hiring  an  additional  post-doctoral  research  associate  (Dr.  Hu  Qi)  to 


assist  in  this  research  effort  were  also  incurred  during  the  year.  Two  candidates  were 
considered  for  this  position,  and  expenses  associated  with  interviews  were  partially 
defrayed  from  this  grant. 

5.  Plan  for  Next  Year:  1  Oct.  1992  -  30  Sept.  1993 

During  the  next  funding  year  of  this  research  effort,  the  cloud  cover  data  and  improved 
parameterizations  of  cloud  coverage  developed  during  the  first  phase  will  be  incorporated 
into  a  mesoscale  meteorology  model.  Model  forecasts  which  utilize  the  observed  cloud 
coverage  and  depth  should  be  improved  relative  to  forecasts  which  crudely  specify  cloud 
properties.  During  this  period,  we  anticipate  running  improved  versions  of  the  MM4 
meteorological  model  on  USAF  supercomputers.  We  plan  to  continue  to  expand  our 
inter  comparison  dataset  on  which  we  analyze  cloud  cover  and  related  meteorology.  In 
addition  to  the  springtime  observation  period,  we  will  be  simulating  a  summertime  period 
in  conjunction  with  RTNEPH  data. 

On  12  October  1992,  Drs.  Walcek  and  Zhang  traveled  to  Phillips  Laboratory  near  Boston 
to  initiate  a  collaboration  with  Air  Force  scientists  involved  with  similar  research  there. 
After  presenting  results  to  the  cloud-analysis  research  group  there.  Dr.  G.  Modica 
familiarized  us  in  the  methods  of  running  the  mesoscale  forecast  model  on  USAF  CRAY- 
2  computer.  Considerable  modifications  to  the  run  procedure  are  required  to  transfer  the 
mesoscale  meteorology  model  from  the  CRAY-YMP  to  the  Phillips  Laboratory 
Supercomputer  Center.  Fortunately,  USAF  scientists  have  experience  with  this  transfer, 
and  can  thus  save  us  considerable  effort  as  our  research  progresses. 


Conference  presentations,  abstracts,  and  journal  manuscripts  resulting  from  this 
research  effort; 

Walcek,  C.  J.,  1992:  Cloud  cover  and  its  reladonship  with  relative  humidity  during  a 
springtime  midlatitude  cyclone:  some  implications  for  climate  models.  Proceedings,  11th 
International  Conference  on  Clouds  and  Precipitation.  Montreal,  Canada,  17-21  August 
1992.  Elsevier  Publishers,  1 1 28- 1131. 

Walcek,  C.  J.,  1992:  Extrapolating  cloud-scale  microphysical,  dynamic,  and  radiative 
processes  to  global  and  climatic  scales:  How  accurately  do  we  know  the  fractional  area 
of  cloud  coverage?.  Workshop  proceedings  of  the  WMO  Cloud  Microphysics  and 
Applications  to  Global  Climate  Change  Workshop.  Toronto,  Canada,  10-14  August 
1992.  In  press. 

Walcek,  C.  J.,  1992:  Cloud  cover  and  its  relationship  with  meteorological  factors  during 
a  springtime  midlatitude  cyclone.  Monthly  Weather  Review,  manuscript  submitted  and 
under  review,  July  1992. 


APPENDIX: 


Manuscript  submitted  and  reviewed  by  the  journal 
MONTHLY  WEATHER  REVIEW 


This  manuscript  has  been  accepted  for  publication  after  minor  revisions. 
Following  some  minor  revisions  (in  progress). 

The  article  will  be  resubmitted  during  the  month  of  November  1992. 
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ABSTRACT 


In  ihis  study.  v\e  compare  vertical  distnhutions  ot  fractional  cloud  coverage  within  -i.'IO 
Nini-  areas  with  co-iivateci  standard  meteorologicai  obserc  ations  over  the  northeast  C.  S.  during  a 
■'pringtirne  midlatitude  cyclone.  Cloud  cover  ('Oser\’ations  are  derived  from  the  I'.  S.  Air  Force 
.'DNF.PH  analysis  (9'  satellite  imagery  and  -'..rr'ace-hased  reports.  Standard  meteorological 
parameters  are  interpolated  from  measurements  using  a  rodrostatic  mesoscale  meteorology  model. 

L'nder  conditions  when  no  camuius  convection  is  possible,  we  find  moderate  to  good 
positive  correlations  between  clouu  cover  and  reiative  humidity  and  large-.scale  venical  veiocii\ . 
.Kid  e\  idence  for  weak  negative  correlations  between  cloud  cover  and  w  md  shear  anu  temperature 
'.apse  rate.  We  pre'-ent  contours  o:  the  mean  fractional  cloud  coverage  observed  at  \anous  reiatn  e 
nimidlties  and  pressures.  ,ind  suggest  ..Igorithms  for  estimating  cloud  coverage  from  humidity 
■ind/or  \enical  velocity.  These  comparisons  suggest  that  cloud  cover  decrea.ses  exponentially  as 
humidity  falls  below  100''^  Relaioe  to  other  iavers  in  the  troposphere,  the  middle  troposphere 
C..'  -  km  above  the  surface)  contains  the  highest  clouu  amounts  at  the  lowest  relative  humidities, 
with  mean  cloud  amounts  of  near  .pOC  humidity  at  b5('>  mb.  .At  the  same  relative  humidity, 
.ireas  with  large-scale  upward  motions  contain  higher  cloud  amounts  than  areas  where  subsidence 
IS  occurring. 

l'nder  corn eciivei\’  unstable  conditions,  we  see  evidence  for  increased  cloud  am.ounts 
iciatme  to  stable  conditions  at  the  >ame  humidity  in  the  upper  troposphere  at  high  humidities.  In 
the  lower  troposphere,  high  huimdit)'  environments  where  convection  is  possible  contain  lower 
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cloud  amounts  relative  to  stable  conditions  at  the  same  relative  humidity,  which  may  result  from  the 
subsidence  of  dr\'  air  into  the  lower  troposphere  under  convectively  unstable  conditions.  These 
findings  must  be  interpreted  with  caution  'ince  we  see  evidence  tor  excessive  water  vapor 
transport  to  the  upper  troposphere  using  the  moisture  convergence-based  convective  mixing 
alizonthm  within  our  dynamic  meteorolotiicai  mterpolation  algorithm.  We  find  evidence  tor  a  bias 
in  our  area-averaged  humidity  estimates  wnen  interpolating  from  sparse,  point  radiosonde 
measuremients  under  unstable  conditions.  More  accurate  methods  of  inferring  cloud  coverage  from 
relative  numidity  will  require  more  accurate  memods  of  inferring  moisture  distnbutions  over  large 
.ireas.  pantcularly  w  ithin  convectively  unstanic  .  reas. 

Thl^  analvsis  suggests  that  there  in  conv;uerable  uncertainty  in  measuring  cloud  cover  and 
other  meteorological  factors  averaged  over  larce  areas,  .suggesting  that  any  formulations  of  cloud 
cover  fraction  are  highly  uncertain.  vDver  a  w  icte  '-ariety  of  meteorological  environments  associated 
with  a  'Dnngtime  midlatitude  cyclone,  mo^t  cloud  cover  parameterizations  used  within 
meteorological,  climate,  .aid  chemical  mouels  of  the  atmosphere  calculate  cloud  amounts  less  than 
reported  by  the  .’DNfPH  observations,  especially  in  the  middle  troposphere,  where  most 
■ilgorithms  specify  zero  cioud  amount::  at  reiative  humidities  below  .SO  -  HiV/r.  Thus,  current  large- 
scale  climate  simulations  or  atmospheric  ciiemicai  nnxieiing  studies  may  significantly  underestimate 
ihe  effects  of  clouds.  More  importantly,  current  climate  models  probably  cannot  adequately 
estimate  ine  potentiallv  Nignificant  chanues  in  cloud  cover  that  can  result  from  small  changes  in 
relative  humidity. 

I.  Introduction 

Microphvsical  ai\i  dvnamical  prixesses  (Vcumng  within  clouds  can  significantly  innuence 
numerous  '.arger-scalc  co  namic.  radiative,  and  chemical  processes  occurring  in  the  troposphere 
fe.g.  .Arakawa  and  .Schunert,  1974;  Ramanathan  et  al..  19S.s;  Walcek  et  al..  1990).  Since  clouds 
,ire  freuiientiy  present,  i:  is  necessan.'  tor  accurate  models  of  tropospheric  climate  or  chemistry  to 
account  tor  their  etfectv  It  is  well  known  mat  clouds  form  when  the  vapor  pressure  of  water 
exceeiis  me  \apor  pressure  that  woiiki  Iv  s.itur.iteii  with  respect  to  liquid  water  or  ice.  Within  any 
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particular  air  mass,  fluctuations  in  temperature  and/or  water  vapor  concentration  can  lead  to  areas 
where  condensation  (clouds)  can  occur  even  though  the  concentration  of  water  averaged  over  me 
air  mass  may  not  be  saturated  at  the  mean  air  mass  temperature. 

In  large-scale  numerical  models  of  the  atm'"^nhere.  chemical  or  meteorological  properties 
can  only  he  e.xpiicitiy  resolved  over  relatively  large  air  masses,  typically  several  lO's  to  lOO's  or 
kilometers  horizontally  and  -1  ()()()  m  vertically,  it  is  not  unusual  to  observe  fluctuations  in  bcith 
temoerature  and  moisture  within  air  masses  of  this  size  due  to  turbulent  motions,  surface 
innomogeneities,  terrain  and  other  factors.  These  perturbations  induce  cloud  formation  on  a  scale 
that  cannot  be  resolved  by  larger-scale  models  of  the  atmosphere.  Therefore,  regional  or  globai- 
'caie  meteorology  models  resort  to  parameterizaiions  of  cioud-induced  radiative,  dynamic,  and 
chemical  processes  resulting  from  sungnd-scale  clouds.  In  most  parametenzations  of  cloud-scale 
processes,  the  heterogeneous  (or  subgnd-scale)  nature  of  cloudiness  is  approximated  by  assuming 
that  a  f  raction  ( f  i  of  each  grid  area  is  occupied  by  clouds.  This  cloud  fraction  is  used  to  apportion 
.1  Mcnificantly  different  cloud  ''forcing  '  into  a  '  •nd-averaged"  fore;.. .  within  gnd  areas  that  contain 
a  mixture  of  clear  and  cloudy  regions.  For  the  results  presented  here,  cloud  cover  fraction  is 
defined  as  the  fraction  of  a  given  horizontal  plane  in  the  atmosphere  where  condensed  water  is 
V  I'lbly  pre.xent.  This  definition  becomes  somewnai  vague  w  hen  "thin"  clouds  are  present,  a 
condition  that  commonly  occurs  in  the  upper  troposphere.  Under  these  conditions,  it  will  be 
difncuit  to  explicitly  qua:  ty  the  fractional  tmea  of  cloud  coverage. 

Many  studies  of  clouds  and  their  effects  on  tropospheric  processes  have  concluded  that 
c\  en  small  cloud  amounts  can  exert  a  significant  intluence  on  larger-scale  processes.  Under  these 
conditions,  the  net  effect  '^f  clouds  on  any  physical  or  chemical  process  w  ill  be  proportional  to  t . 
the  tractional  area  of  cloud  coverage.  Most  models  of  tropospheric  dynamics  assume  that  tne 
fractional  area  of  cloud  coverage  is  determined  by  the  grid-averaged  relative  humidity.  Figure  1 
shi'ws  the  functional  dependence  of  cloud  coverage  within  a  particular  atmosphenc  layer  from  a 
''Urs  ev  of  fomiulations  currentlv  used  by  various  researchers.  .All  formulations  assume  a  "cntical 
relative  humiditv"  of  between  ."'O  -  'M)'  /  .move  which  partially  cloudv  conditions  can  occur.  Below 
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this  critical  humidity,  all  algorithms  specify  totally  clear  skies.  At  humidities  above  the  critical 
humidity,  cloud  fraction  increases  by  differing  functional  dependencies  to  100%  cloud  cover  at 
100%  humidity. 

Figure  1  shows  considerable  differences  between  alternate  formulations  in  assessing  cloud 
coverage  within  current  meteorology  and  climate  models.  At  80%  humidity,  the  NCAR 
Community  Climate  Model  (Williamson  et  al..  1987)  specifies  95%  cloud  cover  (under  stable 
conditions),  the  British  Meteorological  Office  ciimate  model  (Smith.  1990)  uses  0%  cloud  cover, 
while  other  algorithms  specify  cloud  coverage  between  these  e.xtremes. 

The  differences  and  uncenainty  exhibited  in  Fig.  1  have  motivated  researchers  to  evaluate 
and  improve  methods  for  diagnosing  cloud  cover.  Relationships  between  cloud  cover  and  other 
meteorological  factors  can  be  quantified  and  evaluated  using  two  approaches:  (1)  using 
observations  of  cloud  co'.  and  associated  meteorology  (e.  g.  Slingo.  1980;  Sundqvist  et  al., 
1989):  and  (2)  inferring  cloud  cover  from  fine-resolution  dynamic  models  capable  of  explicitly 
resolving  cloud-scale  dynamics  (Xu  and  Krueger.  1991 ). 

The  primary  problem  with  using  observ'ations  to  infer  relationships  between  cloud  cover 
and  other  meteorological  factors  is  the  difficulty  of  measuring  highly-variable  quantities  of  interest 
over  relatively  large  areas.  .Accurate  measurements  capable  of  vertically  resolving  cloud  coverage 
require  high-density  vertical  soundings,  together  with  satellite  and  relatively  expensive  aircraft 
observations.  Due  to  the  expense  as.sociated  with  such  measurement  programs,  there  are  only 
limited  observational  datasets  under  a  few  meteorological  environments  from  which  to  evaluate  and 
improve  cloud  cover  algorithms. 

When  models  capable  of  resolving  cloud-scale  processes  are  used  to  ascertain  cloud 
coverage,  the  accuracy  and  reliability  of  the  modeling  techniques  are  subject  to  an  independent 
evaluation,  which  ultimately  relies  on  observations.  Sensitivities  of  model  calculations  to  the 
model  dimensionality,  turbulence  and  microphysical  parameterizations.  and  other  factors  must  also 
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In  this  study,  we  compare  standard  meteorological  observations  with  satellite-derived  cloud 
coverage  to  ascertain  the  relationships  between  cloud  coverage  and  other  meteorological  factors. 
We  acknowledge  beforehand  that  the  standard  meteorological  parameters  and  the  cloud  cover 
measurements  used  in  this  analysis  contain  a  considerable  level  of  uncertainty  that  is  difficult  to 
quantify.  Howeve'-.  we  feel  that  the  large  size  of  our  comparison  database  will  allow  us  to 
ascertain  the  true  iClationships  that  may  exist  between  cloud  cover  and  other  meteorological 
factors.  With  a  more  accurate,  observation  e-based  method  of  inferring  cloud  cover,  the  global- 
scale  influence  of  cloud  microphysical  and  dynamical  proces.ses  can  be  more  accurately  treated  in 
larger-scale  studies  of  meteorology  and  atmospheric  chemistn,'. 

2.  Cloud  cover  ob.servations 

Numerous  cloud  cover  datasets  have  been  collected  and  analyzed  over  the  past  several 
decades.  Hughes  ( 19841  provides  a  summary  of  the  characteristics  of  cloud  climatologies  available 
in  the  early  198()'s.  The  earliest  cloud  climatologies  were  compiled  solely  from  surface-derived 
obsen.  ations,  and  were  often  aggregated  in  time  at  various  latitudes  and  time  of  day.  Cloud 
ob.servations  in  these  early  archives  were  often  composed  of  once-per-day  observations  averaged 
over  relatively  long  time  frames  into  seasonal  or  monthly  data  at  a  resolution  of  5  -  10°  latitude  or 
longitude.  Virtually  all  of  the  cloud  cover  data.sets  reviewed  by  Hughes  ( 1984)  w'ere  designed  for 
ti.se  in  textbooks  or  climate  studies  requiring  relatively  coarse  temporal  and  spatial  resolutions. 

For  short-term,  regional-scale  meteorological  and  chemical  studies,  instantaneous 
distributions  of  cloud  conditions  are  required  throughout  a  simulation  period,  updated  as  frequently 
as  the  underlying  meteorology  is  changing.  These  models  cannot  readily  u.se  cloud  cover  estimates 
that  are  temporally  averaged  over  large  areas  encompassing  numerous  meteorological 
environments.  These  models  must  use  cloud  coverage  estimates  based  on  shorter-tenn 
observations. 

The  United  .States  .Air  Force  Environmental  Technical  .Applications  Center  has  been 
receivina  and  storing  Air  Force  Global  Weather  Central  (AFGW'C)  cloud  data  since  January  1971. 
From  1971  to  1983. \FGWC  used  an  operational  real-time  three-dimensional  analysis  of  cloud 
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cover  referred  to  as  3D.\EPH  (Fye.  1978).  The  3DNEPH  is  a  global  analysis  of  cloud  cover  that 
uses  surface-based  and  aircraft  reports,  together  with  visual  and  infrared  satellite  imagery  to 
produce  3-D  cloud  cover  information  on  a  routine  basis.  During  periods  when  satellite  or  surface 
data  are  lacking,  clouds  are  inferred  from  rawinsc....e  temperatures  and  dew'  points.  The  data  are 
ended  onto  a  polar-stereographic  global  grid  with  15  vertical  layers  between  the  surface  and  ~16 
km  above  the  surface.  Horizontally,  the  grid  size  varies  from  ~25  km  near  the  equator  to  -60  km 
at  the  poies.  The  3DNEPH  stores  cloud  cover  information  every  three  hours. 

.Any  cloud  cover  databa.se  will  contain  a  level  of  uncertainty  that  is  difficult  to  explicitly 
evaluate.  The  3DNEPH  assimilates  virtually  all  routinely-available  cloud  observations  into  the 
cloud  archive,  making  an  independent  evaluation  impossible.  Numerous  unevaluated  algorithms 
are  used  to  consolidate  surface-ba.sed  observations  together  with  visible  and  infrared  satellite 
imageiA'.  Hughes  and  Henderson-Sellers  1 1985)  penormed  an  analysis  of  a  later  version  of  the  .Air 
Force  cloud  archive  (RTNEPH)  for  1979.  and  although  numerous  areas  of  obvious  but  minor 
errors  were  discovered,  they  found  that  the  RTNEPH  observations  were  generally  reliable  and  in 
good  agreement  with  known  features  of  tropospheric  meteorology.  Problems  were  found  when 
satellite  data  were  gathered  over  highly  variable  backgrounds  or  backgrounds  with  snow  or  sea  ice. 
.Al.so.  penods  of  missing  data  tu'e  not  uncommon  in  the  data,  although  they  are  identified. 

In  this  study,  we  use  five  noon-time  .spring  periods  analyzed  over  the  northeast  U.  S.  by 
the  3DNEPH.  .At  noon,  we  expect  the  mteximum  utilization  of  aircraft  and  surface  repons  together 
with  visual  and  infrared  satellite  imagery  by  the  3DNEPH  analysis.  Fig.  2  shows  the  cloud 
coverage  within  ~(320  km)-  areas  (box  shown  in  Fig.  2)  in  the  layer  800  -  730  mb  at  local  noon, 
23  .Apnl  1981.  The  (320  km)-  averaging  area  corresponds  to  the  finest  "resolution"  of  the 
corresponding  meteorology  observ'ations  used  in  this  analysis.  During  this  five-day  period,  a 
relatively  intense  midlatitude  cyclone  developed  and  traversed  the  domain  shown  in  Fig.  2.  .At  the 
time  of  Fig.  2.  the  cyclone  was  situated  near  the  center  of  the  domain,  and  the  cloud  cover  shows  a 
warm  frontal  region  o\  er  the  Great  Lakes,  and  a  cold  front  extending  from  Pennsylvania  to  Texas. 
This  period  and  domain  were  chosen  for  analysis  since  there  are  a  wide  variety  of  meteorological 
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conditions  present  from  which  we  can  determine  methods  to  evaluate  cloud  coverage  that  may  be 
generally  applicable. 

3.  Standard  meteorology  observations 

Temperature,  moisture  and  dynamical  data  used  in  this  analysis  are  taken  from  observations 
and  spatially  and  temporally  interpolated  onto  an  (80  km)^  Lambert-conformal  grid  using  a 
hydrostatic  meso'-cale  meteorology  model.  Observations  are  derived  from  the  National 
.Meteorological  Cc.iter  global  meteorological  analysis  and  further  enhanced  using  3-hourly  surface 
observations  and  12-hourly  vertical  rawinsonde  measurements.  Vertically,  these  meteorological 
data  encompass  the  surface  and  100  mb  pressure  surface  (~16  km),  and  horizontally  they  span  the 
contiguous  United  States.  The  vertical  gnd  size  of  the  meteorology  data  is  ~80  m  near  the 
surface,  and  on  the  order  of  a  kilometer  or  more  aloft.  These  observations  are  provided  as  initial 
and  boundary  conditions  to  the  NCAR  mes^^  :ale  meteorological  model  (MM4  -  Anthes  and 
Warner,  1978).  During  model  e.xecution.  observ'ations  are  incorporated  into  the  model  calculations 
in  regions  near  observation  locations.  Differences  between  observed  and  calculated  temperatures, 
humidities  and  wind  speeds  are  continuously  minimized  through  the  use  of  additional  tendency 
terms  in  the  momentum,  moisture,  and  thermodynamic  equations  which  "nudge"  the  calculation 
towards  the  observations.  The  mesoscale  meteorological  modeling  system  is  not  used  to  "predict" 
or  "forecast"  meteorology  for  this  study.  Rather,  the  model  provides  a  method  to  enhance  the 
spatial  and  temporal  resolution  of  obseiv'ed  meteorology  by  interpolating  the  observations  in  a 
dynamically  and  physically  reasonably  manner.  Thus,  model  calculations  agree  closely  with 
ob.servations  when  and  where  observations  are  available,  and  w  hen  no  observations  are  available, 
the  meteorological  data  are  dynamically  consistent.  .Meteorology  interpolated  and  analyzed  from 
observations  in  this  manner  provides  data  of  superior  temporal  and  spatial  re.solution  relative  to 
"raw",  point  observations.  Since  virtually  all  available  measurements  are  used  to  "nudge"  the 
model  calculations,  the  temperatures,  moistures,  and  dynamic  vanables  resulting  from  this  tmalysis 
agree  with  the  measurements  to  w  ithin  the  instrument  uncertainty.  The  uncenainty  in  our  estimates 
of  the  meteorological  state  of  the  atmosphere  in  areas  or  times  w  here  no  observ'ations  are  available 
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is  of  course  unknown.  Stauffer  and  Seaman  (1990)  provide  an  analysis  of  the  effects  of 
observation  assimilation  on  the  quality  of  grided  observational  datasets. 

.Any  method  of  extrapolating  “point  ’  measurements  to  larger  averaging  areas  will  introduce 
an  uncenainty  into  the  area-averaged  "obser\'ations”.  With  current  observations,  it  is  virtually 
impossible  to  quantify  the  magnitude  of  this  "scaling-up”  uncertainty.  For  a  particular  radiosonde 
launch,  if  a  balloon-borne  probe  rises  through  a  cloudy  region  of  a  layer  of  the  atmosphere 
containing  a  mixture  of  clear  and  cloudy  areas,  it  will  report  100%  humidity  for  that  layer  when  the 
mean  relative  humidity  in  the  layer  would  be  well  below  i(X)%.  How  often  these  conditions  occur, 
and  just  how  “noisy"  the  temperature  and  w  ater  vapor  fields  are  within  areas  the  size  of  larger-scale 
model  grid  areas  is  a  highly  variable  and  uncertain  quantity.  .Measurement  programs  to  quantify 
these  uncertainties  would  require  high  time  and  space  resolution  vertical  soundings  of  the 
atmosphere  under  a  wide  variety  of  meteorological  conditions.  Until  such  measurements  are  made, 
there  will  be  significant  and  difficult  to  quantify  uncertainties  in  "measurements”  of  any 
meteorological  quantity  averaged  over  areas  larger  than  the  “footprint”  area  within  which  a 
measurement  device  operates. 

Fig.  3  shows  the  relative  humidity  in  the  8(X)  -  730  mb  layer  interpolated  from  observations 
using  the  mesoscale  meteorology  model  described  above.  The  temperature  and  moisture 
calculations  are  aggregated  into  overlapping  (320  kmi-  areas,  representing  a  4  x  4  average  of  the 
80  km  grid  used  by  the  .VI M4.  The  domain  shown  in  Figs.  2  -  3  represents  approximately  one-half 
of  the  domain  simulated  by  the  mesoscale  meteorological  analysis. 

In  Figs.  2-3.  we  have  mapped  both  the  cloud  cover  and  humidity  data  onto  the  identical 
320  X  320  km  grid  by  area-averaging  the  “raw"  temperature,  moisture  or  cloud  cover  data. 
Vertically,  the  cloud  cover  and  mesoscale  interpolation  model  grids  are  slightly  different,  and  the 
cloud  cover  grid  was  transposed  onto  the  meteorology  model  pressure-based  coordinate  using  a 
cloud  volume  con.serving  mapping.  As  a  result,  a  total  of  5600  data  points  at  15  tropospheric 
levels  are  available  for  comparison  during  the  five  noon-time  periods  studied. 
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4.  Cloud  cover  under  stable  conditions 

We  expect  cloud  coverage  to  be  influenced  by  the  presence  of  buoyancy-induced 
convection  that  is  not  resolved  by  the  coarse  resolution  of  the  meteorological  observations  used  in 
this  analysis.  For  the  first  portion  of  this  analysis,  we  consider  only  areas  where  cumulus 
convection  cannot  occur.  Local  stability  (dT/dz)  at  any  point  in  a  sounding  is  not  a  sufficient 
indicator  of  the  presence  of  convective  activity,  since  convection  can  often  penetrate  into 
atmospheric  layers  that  are  absolutely  stable  with  respect  to  venical  perturbations  (Stull.  1991).  In 
order  to  define  areas  ;>  layers  where  buoyancy-mduced  convection  can  occur,  we  provide  a  1  m 
s  '  “Dush"  to  air  with  a  slightly  higher  temperature  and  moisture  content  from  each  point  on  a 
vertical  sounding.  In  a  conditionally  unstable  environment,  the  “pushed”  parcel  will  accelerate 
upwards.  Ignoring  frictional  forces  and  pressure  perturbations,  the  parcel  velocity  at  levels  above 
the  layer  where  it  is  perturbed  can  be  obtained  by  integrating  the  vertical  equation  of  motion  for  a 
parcel  rising  under  the  influence  of  buoyancy  accelerations: 
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where  u'  is  the  parcel  venical  velocity,  is  the  virtual  temperature  of  the  rising  parcel.  is  the 
virtual  temperature  of  the  sunrounding  environment  through  which  the  parcel  rises,  and  g  is  the 
gravitational  acceleration.  The  condensed  water  content  of  the  parcel  (t// )  is  the  total  water  content 
(assumed  to  remain  constant)  of  the  parcel  minus  the  saturated  vapor  mixing  ratio  at  any  level 
above  the  lifting  condensation  level.  For  this  study,  we  define  a  given  layer  to  have  a  potential  for 
convective  clouds  if  parcels  are  capable  of  rising  to  that  layer  from  any  layer  below  under  the 
influence  of  buoyant  forces.  These  areas  have  been  neglected  in  the  following  analysis  in  an 
attempt  to  ascenain  cloud  cover  under  stable  conditions  only.  We  found  between  1000  and  SOCK) 
stable  grid  areas  at  each  of  the  15  tropospheric  levels  during  the  five  noontime  periods  considered. 
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a.  Cloud  coverage  and  relative  humidity' 

The  most  commonly  used  indicator  of  cloud  coverage  is  th.  mean  relative  humidity  of  an 
air  mass.  Others  (Slingo,  1987)  allow  the  local  potential  temperature  lapse  rate  and  venical 
velocity  to  influence  cloudiness.  Here  we  hypothesize  that  venical  shear  of  the  horizontal  wind 
may  also  affect  cloud  cover.  In  order  to  evaluate  whether  any  of  these  parameters  are  correlated 
with  cloud  coverage,  u  e  calculate  the  correlation  coefficient  for  the  best-fitting  linear  regression 
between  cloud  cover  and  each  of  these  meteorological  factors.  While  we  do  not  expect  that  any  of 
these  parameters  will  be  linearly  related  to  cloud  amount,  we  use  this  correlation  analysis  only  as  a 
means  to  initially  idennfy  how-  strongly  each  of  these  factors  correlates  with  cloudiness. 

Fig.  4  shows  the  vertical  distribution  of  the  correlation  coefficient  for  the  best-fit  linear 
relationship  between  clo  'd  cover  and  humidity,  temperature  lapse  rate,  wind  shear,  and  venical 
\elocity.  Correlation  coefficients  with  magnitudes  greater  than  0.8  are  considered  “excellent" 
according  to  standarc  tatistical  texts,  .vhile  correlations  between  0.6  -  0.8  are  considered 
"good",  and  correlations  between  0.4  -  0.6  tu-e  only  "moderate".  Correlation  coefficients  less  than 
0.4  indicate  poor  or  no  relationship  between  ptu'ameters.  This  figure  shows  that  cloud  cover  is 
most  strongly  correlated  (positively)  with  relative  humidity  over  most  tropospheric  layers, 
followed  by  vertical  \elocity.  which  shows  moderate  to  good  positive  correlations  at  most 
tropospheric  levels.  There  is  also  evidence  for  a  weaker  negative  correlation  between  cloud  cover 
and  both  wind  shear  and  potential  temperaiuie  lapse  rate. 

In  the  upper  troposphere,  we  find  almost  no  correlation  between  cloud  cover  and  other 
meteorological  variables.  This  results  from  two  factors:  ( 1 )  cloud  amount  and  relative  humidity 
are  relatively  low  in  this  layer  (-10-20%  cover  at  10-40%  humidity);  and  (2)  measurement 
uncertainty  is  most  likely  significantly  greater  than  the  mean  humidities  and  cloud  cover.  It  is  well 
known  that  water  vapor  concentrations  are  notoriously  difficult  to  measure  at  heights  above  the  300 
mb  pressure  surface,  and  often  the  standard  radiosonde  humidities  or  dew  points  are  not  even 
reported  at  these  heights.  In  addition,  the  highest  "layers"  in  this  analysis  dataset  are  actually 
averaged  over  3  -  4  km  depths  of  the  atmosphere,  and  therefore  the  "cloud  area  fraction"  defined 
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above  actually  becomes  more  like  a  "cloud  volume  fraction”.  Both  relative  humidity  and  cloud 
fraction  used  in  this  analysis  for  the  layers  above  the  300  mb  pressure  surface  are  highly  uncertain, 
and  any  trends  reported  here  must  be  interpreted  with  caution  until  better  measurements  are 
available. 

In  agreement  with  previous  assumptions  of  factors  that  affect  cloud  coverage.  Fig.  4 
suggests  that  cloud  cover  is  most  strongly  correlated  with  relative  humidity.  Using  the 
observations  shown  in  Fig.  2-3,  we  now  further  investigate  the  relationship  between  cloud  cover 
and  relative  humidity.  Data  shown  in  these  two  figures  together  with  observations  from  4 
additional  days  comprise  several  thousand  (320  kmi-  tu'eas  where  we  have  concurrent  observations 
of  both  cioud  cover  and  relative  humidity.  Fig.  5  shows  the  3DNEPH  cloud  cover  in  the  layer 
between  800-730  mb  plotted  as  a  function  of  the  interpolated  relative  humidity  observations  at  over 
23(X)  areas  where  there  was  no  potentiai  for  convection. 

.A  high  degree  of  scatter  is  immediately  evident  in  this  comparison.  The  correlation 
coefficient  for  the  best  linear  relationship  between  cloud  cover  and  relative  humidity  at  this  level  is 
a  relatively  poor  0.64.  This  considerable  scatter  can  be  explained  through  two  hypotheses;  (1) 
cloud  amount  is  not  strongly  related  to  relative  humidity  alone;  and/or  (2).  our  ability  to  measure 
both  cloud  amount  and  relative  humidity  over  large  areas  is  very  poor.  While  it  is  difficult  to 
explicitly  evaluate  the  uncenainties  in  the  meteorological  ptmimeters  presented  in  this  study  without 
canning  out  an  exhaustive  field  measurement  program,  these  uncertainties  are  most  likely  very 
large,  and  as  a  first  approximation  probably  approach  the  ranges  of  the  scatter  shown  in  Fig.  5.  i. 
e.  20-30*^  uncertainty  in  cloud  cover,  and  a  similar  range  for  uncertainty  in  relative  humidity. 
However,  this  high  uncertainty  does  not  imply  that  existing  measurements  cannot  be  used  to  assess 
relationships  between  cloud  cover  and  related  meteorolo^.cal  factors.  If  the  sample  size  for  this 
study  is  large,  and  covers  a  wide  range  of  meteorological  environments,  then  trends  in  cloud  cover 
may  he  discernible,  and  the  functional  dependence  of  cloud  cover  on  related  meteorological 
parameters  may  be  at  least  qualitatively  revealed. 
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Other  factors  must  be  considered  when  viewing  scatter  diagrams  like  Fig.  5.  Some  of  the 
scatter  shown  in  Fig.  5  is  misleading  since  there  ore  several  hundred  points  with  zero  cloud  amount 
at  humidities  below  about  45%  that  are  literally  “on  top"  of  one  another.  Similarly,  there  are  10  - 
20  coincident  points  at  the  100%  humidity  and  100%  cloud  fraction  that  appear  as  a  single  point  on 
this  scatter  diagram,  with  as  much  significance  to  the  viewer  as  individual  outlying  points 
elsewhere  on  the  figure.  In  addition,  factors  other  than  relative  humidity  may  influence  cloud 
cover,  and  thus  contribute  to  the  .scatter.  However,  a  majority  of  this  scatter  is  most  likely 
attributed  to  the  considerable  uncenainty  in  measuring  both  relative  humidity  and  cloud  cover  over 
the  large  areas  considered  in  this  study. 

The  problem  with  obser\'ing  cloud  fraction  or  relative  humidity  over  large  areas  is 
analogous  to  an  observer  using  a  badly  out-of-focus  telescope  to  find  the  location  of  a  star  in  the 
night  sky.  While  the  telescope  obsers'er  might  see  a  light  fuzzy  image  surrounding  the  true  location 
of  the  object  in  the  sky.  the  obserser  could  not  exactly  pinpoint  the  object's  location.  However, 
the  observer  could  plausibly  hypothesize  that  the  true  location  of  the  point  of  light  is  somewhere 
near  the  middle  of  the  fuzzy  image  revealed  by  the  telescope.  Using  similar  reasoning,  we 
hypothesize  that  if  cloud  coverage  is  related  to  relative  humidity,  then  that  relationship  should  fall 
within  the  range  of  (and  most  likely  near  the  middle  of)  the  .scatter  shown  in  Fig.  5.  In  assessing 
trends  in  highly  uncertain  measurements  shown  in  Fig.  5.  we  aggregate  the  observations  into  5% 
relative  humidity  increments,  and  then  average  the  cloud  coverage  within  these  restricted  humidits' 
ranges.  Using  this  averaging  technique,  trends  become  apparent  in  the  highly  scattered  and 
uncertain  observations.  The  average  and  standard  deviation  of  the  cloud  coverage  w'ithin  5% 
relative  humidity  increments  are  shown  as  a  curve  with  error  bra's  on  Fig.  5.  As  expected,  cloud 
amount  increases  as  humidity  increases.  At  humidities  between  20  -  40%.  there  is  10  -  20%  cloud 
cover  on  average. 

This  process  is  repeated  at  all  tropospheric  levels  to  obtain  the  average  cloud  cover  within 
each  layer  at  any  particular  relative  humidity.  Fig.  6a  shows  the  average  cloud  cover  within  (320 
km)-  stable  areas  as  a  function  of  the  layer  relative  humidity  and  pressure.  .At  a  particular  relative 
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humidity,  cloud  amounts  are  greatest  in  the  800  -  600  mb  layer  of  the  troposphere,  a  trend  that  is 
consistent  with  earlier  approximations  (Buriez  et  al..  1988;  Geleyn  et  al.,  1982).  The  highest 
cloud  amounts  occur  under  high  humidities  at  900  -  800  mb.  but  this  figure  shows  that  10-209^ 
cloua  coverage  occurs  at  humidities  as  low  as  20%,  in  contrast  to  the  formulations  shown  in  Fig.  1 
which  ail  specify  zero  cloud  cover  at  humidities  below  50  -  80%.  Standard  deviations  of  cloud 
coverage  within  these  5%  humidity  increments  (Fig.  6b)  fall  between  20  -  30%  in  absolute  cloud 
cover,  which  in  many  cases  is  greater  than  the  mean  cloud  cover. 

The  trends  in  the  average  cloud  amount  shown  in  Figs.  5  -  6  suggest  that  fractional  area  of 
cloud  coverage  decreases  exponentially  as  relative  humidity  falls  below  100%.  and  that  there  is  no 
clear  'critical  relative  humidity"  where  cloud  coverage  is  always  zero.  The  trends  in  the  average 
cloud  amount  shown  in  Figs.  5  -  6  suggest  the  following  approximation  for  cloud  amount  f  as  a 
function  of  relative  humidity  Rli  {Rh<  \ ): 


f  =  exp 


Rh-l\ 
1  -  Rh,  J 


(2) 


where 


1  -  Rh. 


i  0.2  +  ^,  (T<0.75 

'  1.8(1 -(T).  (j>0.75 


'.3) 


where  dis  the  pressure  relative  to  the  surface  pressure  (P/P?).  The  1 1  -  Rh^)  term  in  Eq.  2  -  3  is 
qualitatively  similar  to  the  “critical  humidity"  used  in  previous  cloud  cover  formulations,  although 
here  it  represents  the  relative  humidity  depression  from  1()()%  where  cloud  amount  decreases  to 
37%  (e  y).  Figure  7  shows  the  value  of  ( 1  -  Rh^)  in  Eq.  (2)  that  yields  the  minimum  root  mean 
square  difference  between  observ'ed  cloud  amount  and  cloud  amount  calculated  using  Eq.  (2). 
Also  shown  on  this  figure  are  contours  of  the  root  mean  square  difference  between  cloud  amount 
calculated  using  Eq.  (2)  and  observ'ed  cloud  amount  as  a  function  of  the  values  chosen  for  the  ( 1  - 
/?/;,.)  parameter.  Individual  data  points  represent  the  value  of  ( 1  -  Rh^)  which  produces  the 
minimum  error,  and  'he  two  linear  segments  plotted  show  Eq.  (3),  a  simple  formula  that  attempts 
to  stav  near  the  areas  ^  f  lowest  error  at  each  tropospheric  level.  Using  Eq.  (2)  to  calculate  cloud 
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cover  from  relative  humidity  [both  averaged  over  (320  km)^  areas]  produced  cloud  cover  estimates 
tiiai  on  average  contained  a  root  mean  square  difference  of  7  -  25  percentage  points  from  the 
3DNEPH  obseiA'ations,  depending  on  the  tropospheric  level  considered. 

h.  Cloud  coverage  and  vertical  velocity 

As  shown  in  Fig.  4,  vertical  velocity  is  another  parameter  which  influences  cloud  coverage 
almost  as  significantly  as  relative  humidity.  Some  of  the  scatter  about  the  mean  shown  in  Figs.  5  - 
6  may  result  from  deviations  in  vertical  velocity  within  a  restricted  relative  humidity  range.  In 
order  to  estimate  the  effects  of  vertical  velocity  on  cloud  coverage,  we  further  subdivide 
observations  into  1  cm  s  '  vertical  \  elocity  increments,  and  investigate  relationships  between 
vertical  velocity  and  cloud  cover.  Fig.  8  sr.ows  the  observed  correlation  between  cloud  cover  and 
\  enical  velocity  in  the  800  -  730  mb  layer  with  relative  humidity  restricted  to  50  -  55Cc.  While  the 
correlation  between  venical  velocity  and  cloud  cover  is  weak  at  this  particular  atmospheric  level, 
there  is  a  discernible  trend  in  the  cloud  cover  if  observations  are  aggregated  into  1  cm  S"' 
increments.  In  order  to  facilitate  the  identification  of  trends,  we  have  superimposed  the  nearly 
linear  trend  in  aggregated  cloud  cover  observations  over  the  .scattered  observations.  This  figure 
shows  that  some  of  the  scatter  in  the  50  -  55''f  humidity  range  show  n  in  Fig.  5  can  be  attributed  to 
vanaiions  of  vertical  velocity.  .A  bulk  of  the  low  cloud  amounts  observ'ed  at  this  humidity  have 
venical  velocities  <- 1  cm  s  ',  while  there  are  very'  few  clear  skies  (  /  =  0)  when  the  vertical  velocity 
is  positive. 

In  order  to  quantify  the  influence  of  vertical  velocity  on  cloud  cover,  we  calculate  the 
■'critical  humidity"  term  (1  -  Rhc  in  Eq.  2  that  yields  the  lowest  root  mean  square  differences 
between  me  observ'ed  cloud  cover  and  cloud  cover  calculated  using  Eq.  2  using  observations 
confined  to  1  cm  s  '  vertical  velocity  increments.  Fig.  9  shows  the  functional  dependence  of  this 
critical  humidity  depression  as  a  function  of  vertical  velocity  for  the  tw'o  tropospheric  layers. 
These  data  show  an  approximately  linear  relationship  that  is  reproduced  at  all  other  tropospheric 
levels.  The  fractional  area  of  cloud  coveraue  when  the  vertical  velocitv  is  close  to  zero  is  best 
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approximated  by  Eq.  3.  and  Fig.  lOa  shows  the  vertical  variation  oi  ihe  slope  of  the  best-fitting 
line  through  the  data  shown  in  Fig.  9  for  all  tropospheric  levels.  Therefore,  if  both  relative 
humidity  {Rh\  and  vertical  velocity  uv  -  cm  s  ‘  \  are  used  to  estimate  cloud  coverage,  we  propose 
the  following  formulations: 

I  1  I 

'■=expt - f - : - 7>  ,  (4) 

1  maxf(l-/?/t  )  +  ().lcTMv.O,()()l]j 

Again,  cris  the  pressure  relative  to  the  surface  pressure  (P/Ps),  and  the  ( 1  -  Rh,,)  term  in  the 
denominator  is  specified  in  Eq.  3.  .-\s  shown  in  Fig.  idb.  using  both  humidity  and  vertical 
velocity  lEqs.  3-4)  to  specify  cloud  cover  reduces  the  root  mean  square  deviations  from  the 
3DNEPH  cloud  cover  by  less  than  10  percent  relative  to  using  only  reiaiive  humidity  (Eq.  2  -  3). 
The  calculations  of  Eq.  4  annlicd  to  one  humidity  range  and  level  are  al.so  plotted  as  the  curve 
identified  in  Fig.  8.  For  th  ^articular  level  and  relative  humidity  range,  the  parameterization 
produces  a  slight  overestimate  of  cloud  cover,  although  the  sensitivity  of  the  cloud  cover  to 
changes  :n  venical  velocity  tu'C  well  simulated.  .At  other  relative  humidity  intervals  and  levels.  Eq. 
4  produces  better  or  worse  agreement  with  the  trends  revealed  by  these  observations,  and  on 
average,  the  formulation  suggested  above  produces  the  best  overall  approximation  lore 
complicated  fonnuiations  employing  additional  free  parameters  were  investigated,  but  the 
additional  complexity  of  the  formulations  did  not  significantly  reduce  the  differences  between 
observed  and  calculated  cloud  cover  relative  to  the  formulations  described  above. 

5.  Cloud  cover  under  convectively  un.stable  conditions 

We  now  consider  regions  where  buoyancy-induced  convection  can  intluence  cloud 
coverage.  .As  noted  above,  convectively  unstable  areas  are  defined  as  regions  in  the  atmosphere 
w  here  p.ircels  that  are  slightly  less  dense  than  their  surrounding  environment  can  rise  to  after  being 
given  a  minor  impulse.  Fig.  1 1  shows  the  average  cloud  cover  within  (320  kmi-  unstable  areas  as 
a  function  of  the  layer  relative  humidity  and  pressure.  This  figure  is  qualitatively  similar  to  Fig.  6 
showine  cloud  cover  under  stable  conditions;  i.  e.  cloud  amount  maximizes  in  the  middle 
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troposphere.  Fig.  1  Ih  shows  the  difference  between  the  stable  and  unstable  cloud  covers  /liable ' 
'  ansiabio  Same  relative  humidity  and  pressure.  Since  the  observations  from  which  these 

figures  are  constructed  contain  a  high  level  of  scatter,  some  of  the  differences  between  stable  and 
unstable  cloud  covers  may  not  be  statistically  significant.  Therefore,  we  perform  a  chi-.square 
analysis  of  these  differences.  Relative  numiditv  and  cloud  cover  data  pairs  were  aggregated  into 
cloud  cover  increments  at  each  b'".'  b.umiditv  increment  and  pressure  level.  The  unstable  and 
'table  cloud  cover  frequency  distribuiioi':'  were  then  bjected  to  a  chi-square  confidence  limits 
test.  Differences  between  stable  and  unstable  cioud  cover  that  are  significant  at  greater  than  the 
confidence  limit  are  denoted  by  neaw  \  >  cn  Fig.  1  lb.  and  differences  that  are  significant  at 
treater  than  the  limit  are  denoted  b\  nain  \  Differences  that  are  not  close  to  an  "x"  on  Fig. 
.  lb  are  usually  small  and  may  not  he  statisticaio  'ignificant. 

Fig.  II  shows  [hat  under  relat:\eiy  -...mid  convectively  unstable  conditions,  we  see 
evidence  for  increa.scd  cloud  amounts  re;at;\e  i  able  conditions  at  the  same  humidity  in  the  upper 
troposphere.  In  the  lower  troposphere.  :'.:gn  liumidity  environments  where  convection  is  possible 
contain  lower  cloud  amounts  relative  to  sranle  Cvsnditions  at  the  same  relative  humidity,  which  may 
result  from  the  subsidence  of  dry  air  into  the  .ower  troposphere  under  convectively  unstable 
jonditions.  These  differences  seem  nr.y'ically  reasonable.  However,  under  unstable  conditions 
there  is  a  region  of  significantly  lower  cloud  amounts  when  humidities  are  below  65*^  in  the  upper 
troposphere  that  cannot  be  easily  exoiained.  T  his  anomalous  behavior  prompted  a  further 
investigation  of  the  cloud  cover  and  humidity  viata.  Further  analysis  of  the  MM4  interpolation 
system  showed  that  under  unstable  conditions,  the  upper  troposphere  is  calculated  to  become 
rapidly  saturated  when  obsersuitions  are  niH  assimilated  into  the  model  forecast  equations  (i.  e.  the 
MM4  is  run  in  a  'forecast"  mode  aw  a  .  from  the  model  boundaries).  The  MM4  calculates 
humidities  approaching  SO  -  lOO^r  in  me  iayer  rrom  .S()0  -  100  mb  over  the  region  of  convective 
storms  asstviated  with  transient  cold  aiu:  warm  trouts.  Simultaneous  observations  showed  much 
smaller  regions  of  humidities  in  excess  tn  'O',  txttween  pOO  -  .'(K)  mb.  and  above  250  mb.  relative 
iiumidities  never  exceeded  20'T.  As  note;:  prcMotisly,  these  upper-tropospheric  relative  humidities 
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are  notonousiy  unreliable.  However,  the  differences  between  the  calculated  humidities  and 
“observed"  humidities  are  considerable  in  the  regions  of  frontal  convection.  For  the  version  of  the 
model  usea  in  this  analysis,  any  bia.ses  introduced  by  the  model  physics  into  the  calculation  of 
■vater  vapor  that  is  partially  deleted  since  water  vapor  calculations  are  continuously  “nudged" 
towards  the  relatively  dry  observaitions.  De.spite  this  correction,  there  is  a  possibility  that  under 
convectiveiy  unstable  conditions,  there  may  be  a  bias  in  the  .\lM4-produced  moisture  calculations 
that  employ  tne  moisture-convergence  based  convective  parameterization.  If  this  is  the  case,  then 
the  relative  numidiiies  used  to  construct  Fig.  1 1  may  be  higher  than  were  actually  present, 
especially  inder  the  unstable  conditions  when  the  M.VId  convective  parameterization  was 
calcuiatine  r.mher  moistures  than  observed  in  the  upper  troposphere. 

Thus,  we  see  an  unconsistency  between  the  .\1M4  humidity  calculations  and  the  3DNEPH 
cloud  cover  under  unstable  conditions.  Cloud  amounts  typical  of  relatively  low-  humidities  under 
Ntable  conditions  are  being  observ'ed  at  the  higher  humidities  calculated  by  the  MM4  in  unstable 
areas.  Since  observed  humidities  above  about  .M)()  mb  are  highly  uncertain,  and  possibly 
Nignificantly  lower  than  the  .\IM4  calculations  when  it  is  exerci.sed  in  a  predictive  mode,  it  appears 
that  humidity  calculated  using  .V1M4  under  convectiveiy  unstable  conditions  may  be  overestimated. 
In  short,  there  .  .e  very  few  clouds  observed  w  here  the  .\1M4  is  calculating  SO  -  lOO^f  humidities, 
suggesting  that  the  humidities  may  be  too  high. 

Because  of  these  apparent  discrepancies,  we  feel  that  we  cannot  reliably  use  the 
observation-interpolation  system  de.scribed  above  to  derive  relationships  between  cloud  cover  and 
humidity  under  unstable  conditions.  More  accurate  methods  of  correlating  cloud  cover  and  relative 
humidity  or  other  meteorological  factors  w  ill  require  more  viccurate  methods  of  inferring  moisture 
distributions  over  large  .ireas.  particularly  within  convectiveiy  unstable  areas.  This  apparent 
discrepancy  may  also  intluence  humidity  “measurements"  reported  here  under  stable  conditions  in 
the  upper  tronosphere  since  this  is  an  integrated  observation  assimilation  system,  although  these 
errors  would  be  contmuousiy  minimized  in  areas  where  observations  are  available.  However,  as 
noted  earlier,  in  the  upper  tropospiiere  where  even  direct  measurements  of  water  vapor  are  highly 
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uncenain.  any  formulations  representing  relationships  between  cloud  cover  and  humidity  derived 
from  those  measurements  are  open  to  considerable  uncenainty. 

6.  Discussion  and  Conclusions 

In  this  study,  we  have  compared  satellite  observations  of  fractional  cloud  coverage  within 
~<320  km)^  areas  with  related  meteorological  observations  over  the  northeast  U.  S.  during  a 
springtime  midlatitude  cyclone.  Cloud  cover  observations  were  derived  from  the  U.  S.  Air  Force 
3DNEPH  analysis  of  satellite  imagery,  aircraft  repons,  and  surface-based  observations.  Other 
meteorological  measurements  were  interpolated  from  radiosonde  observations  using  a  hydrostatic 
mesoscale  meteorology  model.  Co-located  comparisons  of  the  cloud  cover  with  other 
meteorological  measurements  show  considerable  uncertainty,  although  we  find  moderate 
correlations  between  cloud  cover  and  relative  humidity  and  vertical  velocity.  .An  analysis  of  the 
relationship  between  observed  cloud  cover  and  humidity  suggests  that  cloud  cover  decreases 
exponentially  as  humidity  falls  below  100%.  Relative  to  other  layers  in  the  troposphere,  the 
middle  trooosphere  (700-500  mb)  contains  highc'  :loud  amounts  at  lower  humidities,  with  mean 
cloud  an'.uunts  of  ~30%  near  50%  humidity.  At  the  same  relative  humidity,  cloud  cover  increases 
with  increasing  vertical  velocity,  and  cloud  cover  also  exhibits  weak  negative  correlations  with 
potential  temperature  lapse  rate  and  venical  shear  of  the  horizontal  wind. 

The  analyses  described  above  were  repeated  at  several  different  horizontal  averaging 
resolutions.  The  fundamental  observations  were  horizontally  averaged  from  areas  representing  the 
horizontal  grid  size  of  the  meteorological  interpolation  model  (80  km)  up  to  8  x  8  averages  (640 
km)  of  the  finest  resolution  observations.  There  were  only  minor  quantitative  fluctuations  of  the 
correlations  described  above,  and  the  only  parameters  which  changed  significantly  were  the  root 
mean  square  deviations  and  the  correlation  coefficients  between  formulations  of  cloud  cover  based 
on  other  meteorological  factors.  As  the  horizontal  averaging  area  became  larger,  the  “scatter”  (e.  g. 
see  Fig.  5)  became  smaller  as  numerous  grid  areas  with  deviations  became  merged  into  a  single 
area.  Thus,  the  above-described  methods  for  estimating  cloud  cover  within  a  particular 
atmospheric  layer  should  be  applicable  for  all  model  grid  sizes  >  80  km. 
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onder  conditions  when  buoyancy-driven  vertical  motions  are  possible,  we  find  some 
discrepancies  between  model  calculations  of  relative  humidity  and  the  cloud  cover  observations, 
suggesting  a  bias  in  the  method  of  interpolating  relative  humidity  in  the  upper  troposphere.  In 
particular,  we  .see  low  cloud  amounts  that  are  typical  of  much  drier  conditions  occurring  at  higher 
humidities  under  unstable  conditions,  which  we  attribute  to  excessive  vertical  transport  of  water 
vapor  under  convective  conditions.  We  are  thus  using  the  3DNEPH  cloud  observations  as  an 
evaluation  tool.  In  a  manner  similar  to  Nehrkorn  and  Hoffman  ( 1990),  our  analysis  suggests  that 
cloud  coverage  may  be  a  useful  tool  to  infer  relativ-  imidity,  contrary  to  the  underlying  purpose 
of  most  cloud  cover  algorithms:  to  infer  cloud  coverage  from  other  "easy  to  estimate" 
meteorological  factors.  For  inferring  upper  troposphenc  water  vapor  concentrations,  perhaps 
satellite-derived  cloud  coverage  is  a  parameter  that  is  easier  to  estimate  with  greater  accuracy  than 
relative  humidity. 

Despite  the  high  level  of  uncertainty  present  in  this  analysis  of  cloud  cover,  it  is  obvious 
that  most  parametenzations  of  cloud  coverage  used  by  current  climate  models  calculate  smaller 
cloud  amounts  than  reported  by  the  3DNEPH  observations,  especially  in  middle  tropospheric 
levels  at  relative  humidities  less  than  70  -  80%.  It  is  apparent  that  many  cloud  cover  formulations 
show  n  in  Fig.  1  assume  that  cloud  cover  fractions  are  considerably  lower  than  the  mean  minus  one 
or  more  standiu'd  deviations  from  the  measurements  presented  in  Fig.  5.  The  differences  betw'een 
current  formulations  and  these  observations  is  greatest  at  low  humidities.  Many  of  the 
formulations  for  cloud  cover  used  by  climate  models  are  not  actually  based  on  short-term 
observations  of  cloud  cover  and  relative  humidity.  Rather,  these  formulations  along  with  their 
"critical  humidities"  ,ire  "tuned"  w'ithin  the  context  of  their  host  model  physics  and  dynamics  to 
yield  reasonable  estimates  of  the  global,  long-term  planetary  albedo.  Thus,  cloud  cover 
formulations  are  "backed-out"  of  their  host  model,  and  represent  a  functional  form  which  when 
used  in  the  panicular  climate  model  for  which  they  were  tuned,  yield  a  reasonable  planetary  albedo. 
The  danger  of  this  method  of  assessing  cloud  cover  is  that  other  model  uncertainties  and  even 
errors  can  be  compensated  for  by  an  incorrect  cloud  cover  fonnulation.  For  example,  if  a  model 
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uses  an  optical  depth  for  clouds  that  is  too  great,  then  this  error  can  be  compensated  for  by  making 
the  clouds  smaller  (i.  e.  a  smaller  cloud  fractional  coverage).  Most  cloud  cover  formulations 
appear  to  calculate  cloud  amounts  that  cure  considerably  smaller  than  reported  in  the  3DNEPH 
tu'chive.  yet  these  same  climate  models  claim  to  reproduce  the  global  planetary  albedo  accurately. 
This  suggests  that  other  components  of  climate  models  may  be  uncertain  or  in  need  of 
improvement. 

These  results  suggest  that  current  methods  of  calculating  cloud  coverage  within  large-scale 
climate  simulations  or  atmospheric  chemical  modeling  studies  are  significantly  underestimating 
cloud  amount  at  most  relative  humidities.  .More  importantly,  current  climate  models  probably 
cannot  adequately  estimate  the  potentially  significant  changes  in  cloud  cover  that  can  "-esult  from 
small  changes  in  relative  humidity.  .As  is  clearly  shown  in  Fig.  1.  a  small  change  in  relative 
humidity  can  result  in  large  or  small  changes  in  cloud  coverage,  depending  on  which  cloud  cover 
algorithms  are  used.  This  analysis  of  the  3DNEPH  cloud  archive  shows  an  exponentially 
decreasing  cloud  cover  as  relative  humidity  decreases,  and  no  clear  "critical  relative  humidity’’ 
below  which  there  are  no  clouds.  More  importantly,  according  to  the  trends  revealed  by  this 
analysis,  cloud  amount  probably  changes  as  relative  humidity  changes  at  any  relative  humidity, 
unlike  current  fomtulations.  which  only  allow  cloud  cover  to  change  when  relative  humidity  is 
greater  than  60  -  90‘T. 
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Relative  humidity  {%) 

FIG.  1.  Fractional  cloud  coverage  as  a  function  of  relative  humidity  at  800  mb  according  to 
various  formulations  used  by  meso-  and  global-scale  atmospheric  models. 
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FIG.  4.  Correlation  coefficient  at  various  atmospheric  levels  for  the  best-fit  linear  relationship 
between  cloud  cover  and  relative  humidity,  potential  temperature  lapse  rate,  wind  shear,  and  (320 
km)-  averaged  venical  velocity. 
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FIG.  5.  Fractional  cloud  coverage  as  a  function  of  relative  humidity  at  800-730  mb.  Each 
point  represents  one  (320  kmi^  area  in  the  domain  shown  in  Fig.  3  during  20-24  April  1981.  Only 
gnd  areas  where  no  buoyancy-driven  convection  can  occur  are  considered.  Lines  show  the  mean 
and  standard  deviation  ot  cloud  cover  within  5%  increments  of  relative  humidity. 
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FIG.  6.  (a)  Fractional  cloud  coverage  and  (b)  standard  deviation  ot  fractional  cloud  coverage 
as  a  function  of  relative  humidity  and  pressure  during  20  -  24  April  1981  over  the  northeast  L.  S. 
shown  in  Fitz.  3.  Onlv  areas  where  no  convection  occurs  are  used  to  construct  this  figure. 
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(1-Rhe)  •  relative  humidity  depression  where  /  =  1/e 

FIG.  7.  Contours  of  the  root  mean  square  error  (%  cloud  cover  )  in  calculating  cloud  amount 
from  relative  humidity  using  Eq.  (2)  as  a  function  of  values  for  the  crinca!  humiditv  depression  ( 1- 
/?/iei  and  pressure  level  in  the  atmosphere.  Points  represent  the  values  for  [l-Rfic)  that  yield  the 
minimum  difference  between  observed  and  calculated  cloud  faction.  Solid  curve  shows  Eq.  t3'). 
the  suggested  vertical  variation  of  (\-RhJ. 
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FIG.  8.  Fractional  cloud  coverage  as  a  function  of  vertical  velocity  at  800-730  mb.  relative 
humidity  between  50  -  55  Each  point  represents  one  (320  km)*  area  in  the  domain  shown  in 
Fig.  3  dunng  20-24  .Apnl  1981.  Only  gnd  areas  where  no  buoyancy-driven  convection  can  occur 
are  consiaered.  Line  identified  as  "trend”  shows  the  mean  cloud  cover  within  1  cm  s  ‘  venicai 
velocity  increments.  Line  identified  as  "Eq.  3-4"  shows  cloud  amount  calculated  using  Eq.  3  -  4. 


Vertical  velocity  (cm/s) 

FIG.  y.  Critical  relative  humiditc  depression  that  yields  minimum  root  mean  square  difference 
between  calculated  and  observed  cloud  cover  as  a  function  of  vertical  velocity  at  the  8(K)  -  730  mb 
layer. 
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FIG.  10.  (a)  Vertical  variation  of  the  sensitivity  of  critical  relative  humidity  depression  to 
vertical  velocity  (Rh  change  per  cm  s  '  -  slope  of  best  fitting  line  shown  in  Fig.  9).  Data  points 
plotted  are  the  values  in  Eqs.  4  which  yield  the  minimum  RMS  difference  between  3DNEPH  and 
cloud  cover  calculated  from  relative  humidity  and  vertical  velocity,  (b)  Root  mean  square 
differences  between  cloud  cover  calculated  using  Eqs.  Z  -  4  and  observed  cloud  cover  during  20  - 
25  April  1981. 
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FIG.  11.  (a;  Fractional  cloud  coverage  as  a  function  of  relative  humidity  and  pressure  during 
20  -  24  April  1981  over  the  northeast  L'  S.  (domain  shown  in  Fig.  ?).  Only  areas  where 
buoyancy-driven  convection  can  occur  are  considered  in  this  figure,  (b)  Difference  between  stable 
and  unstable  cloud  cover  (.f stable  '  /unstable-  ‘  F*?-  1  a  function  of  relative  humidity 

and  pressure.  Heavy  x  s  denote  differences  that  are  significant  with  greater  than  95%  confidence, 
lighter  x's  denote  differences  that  are  significant  above  a  90%  confidence  threshold.  .All  other 
areas  contain  differences  that  are  significant  at  less  than  90%  confidence.  Shaded  areas  in  Fig.  1  lb 
show  areas  where  cloud  cover  is  greater  under  convectively  unstable  conditions  relative  to  stable 
conditions  at  the  same  pressure  and  relative  humidity. 


